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We propose a new approximation scheme to solve the Non Perturbative Renormalization Group 
equations and obtain the full momentum dependence of n-point functions. This scheme involves 
an iteration procedure built on an extension of the Local Potential Approximation commonly used 
within the Non Perturbative Renormalization Group. Perturbative and scaling regimes are accu- 
rately reproduced. The method is applied to the calculation of the shift ATc in the transition 
temperature of the weakly repulsive Bose gas, a quantity which is very sensitive to all momenta 
intermediate between these two regions. The leading order result is in agreement with lattice cal- 
culations, albeit with a theoretical uncertainty of about 25%. The next-to-leading order differs by 
about 10% from the best accepted result. 
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In this paper we present a new and systematic method 
of approximation for the Non Perturbative Renormali- 
sation Group (NPRG) (3, 0, , which allows one to get, 
in a simple way, the full momentum dependence of the 
n-point functions. During the last years, the NPRG has 
been applied successfully to a variety of physical prob- 
lems (for reviews, see e.g. 0, |1|)- In most cases, the 
solution of the NPRG equations involves a derivative ex- 
pansion which only allows for the determination of the 
n-point functions and their derivatives at vanishing ex- 
ternal momenta. There are however situations where the 
knowledge of the full momentum dependence of the n- 
point functions is needed. A simple example is that of 
the calculation of the shift ATc of the transition temper- 
ature of a weakly interacting Bose gas 0. This will be 
used here as a test of our new approximation scheme. 

In order to present our method, we consider a scalar 
(j)^ theory in d dimension with 0{N) symmetry: 



The role of AS'^ is to suppress the fluctuations with mo- 
menta p ^ K, while leaving unaffected the modes with 
p ^ K. Thus, typically Rk{p) — > when p <^ k, and 
Rk{p) ~* when p > k. There is a large freedom in the 
choice of Rk{p), abundantly discussed in the literature 
0, 0, S lioll. We have used a cut-off function proposed 
Q [13: 



by Litim [lOj: Rk{p) oc (k^ —p^) —p^), which allows 
many calculations to be done analytically. 

The NPRG equations can be written as an infinite hi- 
erarchy of flow equations describing how the n--point 
functions evolve with the infrared cut-off k P, |2j 0]- 
Here, we write explicitly only the equation for the self- 
energy E and that for the 4-point function F'^'*), which is 
enough for the purpose of explaining our approximation 
scheme. For vanishing fields, {(f) = 0, the flow equation 
for I]jj(K;p) = 5ijYi{K;p) reads 



5ijdKX'(i<,;p) 
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+ ^ W^^)] I d''^; . (1) and that for P^^^].; (k;pi,P2,P3,P4) is given by 



The field (t>{x) has N real components (piix), with i = 
I,--- ,N. 

The starting point of the NPRG is a modification of 
the classical action to which is added 



{ 



(27r)" 



-Pi 



d^R4q)G^{K;q) 



^SM = \j ^R.{p)Mp)M~p)- 



(2) 



- P2)G(K;q+pi +P2) 
xrfctim('^;P3>P4, q~P3~ Pi, ~q) + permutations 

i(K;pi,P2,P3,P4,g, . (4) 



r 



(6) 

ijklmm V 



•Electronic address 
t Electronic address: 
■f Electronic address: 



I blaizot Qspht . saclay. c ea. t'r | 
mendezg@fing.edu. uyj 
nicws@fing.edu. uyj 



More generally, the equations for F^"^ involves all the 
pC™) with m < n < n+2. In these equations, G^^{k;p) — 
p^ + Rk{p) + E(/t,p). These equations are to be solved 
with the boundary condition that as k — > A, with A 
the microscopic scale, the n-point functions take their 
classical values, read on the action (^). In particular. 
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^ijkii^'^Pi'P2,P3,P4) = {u/3){5ijdki + SikSji + SiiSjk)- As 
K — > the n— point functions go to their physical values. 

Clearly, in general, the NPRG hierarchy of equations 
can only be solved approximately. But a virtue of the 
NPRG is precisely to suggest approximations which are 
not easily formulated in other, more conventional, ap- 
proaches. A popular one is the derivative expansion 
[L llj ■ In lowest order it consists in ignoring the momen- 
tum dependence of vertices in the right hand side of the 
flow equation, as well as field renormalization. In this 
approximation, usually referred to as the Local Poten- 
tial Approximation (LPA), the hierarchy collapses into 
a closed equation for the effective potential Vk((/)^). An 
interesting improvement of the LPA, which we refer to as 
the LPA', takes into account a running field renormalisa- 
tion constant and allows for a non trivial anomalous 
dimension, determined from the cut-off dependence of 
Zi^, Tj = — InZ^ IJ. The solution of the LPA' is well 
documented in the literature (see e.g. |5j,i^J). It will be 
used as an input in our approximation scheme. Let us 
emphasize that for massless theories, the derivative ex- 
pansion is an expansion in terms oip^ /k^ and therefore, 
in the physical limit k — » 0, it makes sense only for van- 
ishing external momenta, p — Q. Higher orders in the 
derivative expansion ^] may improve the accuracy of the 
LPA but, again, this concerns only the n-point functions 
or their derivatives at zero momenta. 

In order to get the full momentum dependence of the 
n-point functions, we propose here to solve the hierarchy 
of the NPRG equations through an iteration procedure. 
The iteration of the NPRG equations starting with the 
classical values of the n-point functions as initial input 
reconstructs the usual loop expansion — see e.g. 12. 13i]. 
Thus one expects generically the iterations of the exact 
flow equations to correctly account for the high momen- 
tum behaviour. This is not so however for the low mo- 
mentum region. There, our approximation scheme will 
exploit the fact that the LPA' is a good approximation 
in the limit of vanishing momenta. 

The iteration procedure starts with an initial guess for 
the n-point functions to be used in the right hand side 
of the flow equations. Integrating the flow equation of a 
given n-point function gives then the leading order (LO) 
estimate for that n-point function. Inserting the leading 
order of the n-point functions thus obtained in the right 
hand side of the flow equations and integrating gives then 
the next-to-leading order (NLO) estimate of the n-point 
bmctions. And so on. 

There is no small parameter controling the convergence 
of the process, and the terminology LO, NLO, etc. that 
we just used refers merely to the number of iterations 
involved in the calculation of the n-point function con- 
sidered. Obviously, the calculations become increasingly 
complicated as the number of iterations increases, and it 
is essential that the initial guess be as close as possible 
to the exact solution so that only one or two iterations 
suffice to get an accurate result. The main effort focuses 
then on the construction of such a good initial guess for 



the solution, to which we now turn. 

We shall be guided by a crucial property of the NPRG, 
that of decoupling of the various momentum scales: for 
given external momenta p ^ a, the flow of the n-point 
functions is dominated by internal momenta q ^ k (the 
derivative d^RniQ) limits the range of integration in the 
flow equations to g ^ k), and when k becomes lower 
than p the flow essentially stops (the external momenta 
playing the role of infrared regulators). 

Consider first the 2-point function. A reasonable ini- 
tial guess for the propagator to be used in the r.h.s. of 
the flow equations when calculating the leading order es- 
timate of the n-point functions is the LPA' propagator 

^LPA' i'^'i P) ~ '^K.P^ + ™K + (p) ' where rn^ is the run- 
ning LPA' mass: it is a good approximation when p k, 
and for p ^ k, G~^{k\p) goes quickly to zero and its 
precise form does not matter. 

A good initial ansatz for the other n-point functions 
(with n > 2) will be obtained by solving the flow equa- 
tions for the r*^"-''s, with the following three approxima- 
tions. 

Our first approximation (Al) assumes that for q ^ k, 
and any set of external momenta {pi,P2, ■■■,Pn} we have 

r'")(«:;pi, ...,pn-i + g,p„ ~ q) ~ r<")(«:;pi, ^ ^ 

r(")(«:;pi,...,p„) 

(5) 

This is certainly true for momenta {pi,P2j •■•jPn} whose 
norms are much larger than k, if F^") is a smooth function 
of its arguments. Similarly, for vanishing momenta, and 
assuming again that r(") is a smooth function, the con- 
dition above is equivalent to saying that one can expand 
in powers of which leads in zeroth order to the 

LPA, a good approximation. Based on this assumption, 
we shall set g = in all P*^"^ (n > 2), and factor them 
out of the integral in the r.h.s. of the flow equations. 

The second approximation (A2) concerns the propaga- 
tors in the flow equation, for which we make the replace- 
ments: 

G{p + q)^GLPA'{q)eil-^) (6) 

where a is an adjustable parameter. As an illustration of 
the quality of this approximation, we show in Fig. ^ the 
ratio Ik{p)/Ik{0) where 

Up)^ J ^dMq)G'{n;q)G{K;p + q) (7) 

is the integral wich remains in Eq. Q after approxima- 
tion Al. This ratio, as a function of p^ / k^, looks indeed 
like a step function, with a weak residual k dependence. 
The approximation A2 amounts to truncate severely the 
high momentum tails of the propagators. This causes a 
slight inaccuracy at high momenta, and a dependence of 
the leading order results on the value of a. One may fix 
a so that the inflexion point of the curve in Fig. is at 
ap = K. One then obtains typically a « .9. One can 
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FIG. 1: The function Jk(p)//k(0) as a function of ln(p/K) 
calculated with the leading order propagator at criticality 

also adjust a so that the integral over k of Ik{p) is iden- 
tical to that of /k(0)6(k^ — a^p^). This yields typically 
a w .6. We regard these two possible choices as extremes 
and adopt the value a = .75 ± .15 for our leading order 
estimate. 

The third approximation (A3) concerns the function 
p("+2) jjj ^jjg equation for P'"^. For this we use an ansatz 
inspired by the expressions of the various n-point func- 
tions in the large TV limit 0|. To be specific, consider 
r^^) in the equation for F^*). The approximation A3, 
combined with Al and A2, leads to the result that the 
contribution of F^^' to the r.h.s. of Eq. Q is proportional 
to that of the other terms, the proportionality coefficient 
i^K being only a function of k. The same proportionality 
also holds in the LPA regime, which allows us to use the 
LPA to determine F,^. 

Approximations A1-A3, when applied to Eq. (0J, yield 
the following equation for our initial guess for F*^**): 

d.r%{K;pi,p2,P3,Pi) = li?\0) (1 - F.) 

X |e(K^ - +P2 f)) r'f-l^„{K-pi,p2,0,~pi -P2) 

xr^'|^^(K;P3,P4, -P3 - P4, 0) + permutationsj . 

(8) 

This equation can be solved analytically in terms of the 
solution of the LPA'^^I- This is done by steps, start- 
ing form the momentum domain a^{pi +p2Y,oP'(pi + 
P3)^,a^(pi + PAf' < K^, where it can be verified that 
the solution is that of the LPA' itself. The solution can 
be written explicitely in the various momentum regions 
defined by the 6-functions occuring in Eq. JSJ. 

The leading order result for the self-energy is obtained 
by using the solution of Eq. © together with Glpa' in 
the r.h.s. of Eq. and integrating numerically over n. 



FIG. 2: The leading order function a{x — p/u) at criticality, 
and the integrand of Eq. 1101 as a function of \n{p/u), for 
u ^ 10^'' for which a{x) is independent of u llflj 

The resulting self-energy at criticality, i.e., for the value 
of r that yields a vanishing physical mass, is shown in 
Fig. 121 for N = 2 and d = 3. It has good behaviour at 
both low and high momenta, independently of the value 
of a: In the scaling regime, one recovers the expected 
result p^ + ^{p) oc p'^~^ , where the anomalous dimension 
is 7y ~ 0.043; this is only slightly larger than values deter- 
mined by the most accurate available methods |l5j | . One 
can show that the value of rj coincides with that ob- 
tained within the LPA', but now it results directly from 
the momentum dependence of the self-energy, rather than 
the cut-off dependence of the field renormalization. At 
large momenta, the logarithmic behaviour given by per- 
turbation theory is reproduced. However, the coefficient 
of the logarithm differs by about 10% from the expected 
one. This problem, cured in NLO, finds its origin in the 
approximation A2 which truncates the high momentum 
tails of the propagators. 

We have applied this approximation scheme to the cal- 
culation of the shift ATc of the transition temperature of 
a weakly interacting Bose gas. It has been shown recently 
that ATc is linear in an^^^ where a is the scattering 
length and n the particle density: 



Here T° is the condensation temperature of the ideal gas, 
given by nXl = C(3/2) with = 2n/mF° (m is the 
boson mass), and ATc = Tc — F^ with Fc the transi- 
tion temperature of the interacting system. As shown in 
Ref. iJ^J, the coefficient c can be related to the change 
of the fluctuation of the field described by the 
action CEJ: c = -2567r3 (^(3/2))"^/^ (A((/.2)/iVu), with 
N = 2. The parameters r and w in Q are then related 
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to the scattering length a and the chemical potential /i 
by: u = 967r2a/A2, and r = -2mTn 19]. 

The best numerical estimates for A(0^), and hence for 
c, are those which have been obtained using the lattice 
technique by two groups, with the results: c ~ 1.32±0.02 
[l3| and c = 1.29 ± 0.05 18]. The availabilty of these re- 
sults has turned the calculation of c into a testing ground 
for other non perturbative methods: ex pan sion in 1/iV 
[Tgl I20I ] , optimized perturbation theory |2l| , resummed 
perturbative calculations to high loop orders Note 
that while the latter methods yield critical exponents 
with several significant digits, they predict c with only 
a 10% accuracy. This illustrates the difficulty of getting 
an accurate determination of c using (semi) analytical 
techniques. 

To understand better the origin of the difficulty, let us 
write A(0^) as the following integral |0] 

A(0f) ^ f dx a{x) 
Nu J 27r2 a;2 + a{x) ' ^ ' 

where cr(x) — u^'^T,{p — xu), with S(p) the self-energy 
at criticality, i.e., I](0) = 0. The integrand of Eq. (^UJ), 
at leading order, is shown in Fig. 2. The momentum 
at the minimum defines the typical scale which separates 
the scaling region from the high momentum region where 
perturbation theory applies. The difficulty in getting a 
precise evaluation of the integral H10|l is that it requires 
an accurate determination of S(p) in a large region of 
momenta including the crossover region between two dif- 
ferent physical regimes 0, . 

In leading order, we obtain c « 1.3 ± .3, in agreement 



with lattice results |23|. The large uncertainty comes 
from the freedom in the choice of a that we have dis- 
cussed above. We have also studied the NLO 0]. This 
involves the leading order for F^^^ which requires the 
initial guess for F*^^^. The latter is obtained by following 
the same strategy as for F^''^ The NLO result for c still 
depends on a, however in such a way that one can fix a 
from a criterium of fastest apparent convergence: in fact 
there is a value of a for which the correction to the LO 
result vanishes. This value is a = .83 leading to c = 1.44. 

Results also exist for other values of N , from lattice 
calcultions 0,0] or perturbation theory to 7-loops [2^ . 
For = 1, we get c = 1.20, to be compared with c = 
1.09±0.09 (lattice) or 1.07±0.10 (7-loops); for iV = 3, we 
get c — 1.62, to be compared with 1.43 ± 0.11 (7-loops); 
for iV = 4, we get c = 1.79, to be compared with 1.60±.10 
(lattice) and 1.54±0.11 (7-loops). Our results thus differ 
by 10 to 15% from the best estimates. In its present 
numerical implementation, our method looses accuracy 
in the large N limit, where we get c « 2.8 instead of the 
exact result c = 2.3 [Tgj . 

In summary, we have proposed a simple method to get 
the momentum dependence of the n-point functions us- 
ing the NPRG. Using well motivated approximations we 
have produced a simple solution which provides an excel- 
lent starting point for a more accurate iterative solution, 
as revealed by its application to the calculation of ATc. 
Although there is no small parameter that controls the 
convergence of the iterative procedure, the NLO calcula- 
tion suggests that the scheme is stable. 
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